Multimode Interference: Identifying Channels and Ridges in Quantum Probability 

Distributions 



- 

X 



Ross C. O'ConnelQ 

Department of Physics 
University of Michigan 
Ann Arbor, MI 48104 

Will LoinaJl 

Department of Physics 
Amherst College 
Amherst, MA 01002 

' The multimode interference technique is a simple way to study the interference patterns found in 

I many quantum probability distributions. We demonstrate that this analysis not only explains the 

, existence of so-called "quantum carpets," but can explain the spatial distribution of channels and 

' ridges in the carpets. With an understanding of the factors that govern these channels and ridges we 

have a limited ability to produce a particular pattern of channels and ridges by carefully choosing 
, the weighting coefficients c„. We also use these results to demonstrate why fractional revivals of 

initial wavepackets are themselves composed of many smaller packets. 



m 

> ■ I. INTRODUCTION 

(N ■ 

' The existence of multiple time scales in many simple, one-dimensional quantum systems is both an important 
characteristic that distinguishes them from their classical counterparts and the source of many interesting phenomena 
in the time-evolution of those systems. Let us begin with a wavefunction of the form 

' oo 

5' vI;(:j,t)^^c„Vn(2:)e-^^"*, (1) 

^ ' n=l 

. . . 

^ I ^ , where ipn is the nth eigenfunction, En the energy of the nth eigenfunction, and the c„ are weighting coefficients. The 
r"| most straightforward way to identify the time scales is to treat £'„ as a function of n and perform a Taylor expansion 
d of the energy around some quantum number n, the number around which the weighting coefficients are centered, 

^3 : 



EnKEn+ E'n {n-fi) + -E'l {n - nf + ... (2) 



We identify the time scales as 



The Tj almost always obey the hierarchy Ti ^ T2 <C Tb-.-, as is demonstrated in Appendix IXI Ti corresponds 
to the classical period of a particle of energy En, so we call it Td- At times near T2 the contribution from the 
E'^ {n — n) terms to the phase of each eigenfunction will be negligible. When the aforementioned hierarchy holds, 
the contribution from the Tj terms for j > 3 will also be negligible and the classical behavior governed by Td will 
dominate, approximately reproducing the t fa behavior of the particle. The return of the wavefunction to its t « 
value is called a quantum revival, and we call T2 the revival time, Tr. The literature on the subject of revivals is 
well-developed, including several general articles 0, QJz^, treatments of fractional revivals 0, 0, 0, H^, and 
articles on possible applications of this phenomenon II, M, 0| • Revivals are also analogous to the Talbot effect in 
classical electrodynamics, as is discussed in [rL 1^ Illll2(1 | . 

We wish to explain the behavior of the wavefunction at times comparable to, but not equal to, Tr. At these times, 
in problems with anharmonic spectrum and for wavefunctions "^{x, t) in which the distribution of c„ has central value 
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n and characteristic width An which obey the hierarchy 1 <C An <C n (a semiclassical case), the spacetime diagram of 
the probabihty distribution exhibits a complex pattern of channels and ridg es - what has come to be called a quantum 
carpet. The problem has been a ppr oached with several techniques [13, HI, EM 122, 123 1 ^'^d finds application in study 
of Bose-Einstein condensates [ol l28l *29^. In this paper, we develop the technique of multimode interference \v\ Il8| . 
In section^ we briefly review the multimode interference technique and introduce characteristic velocities. In section 
mil we obtain, in the context of the WKB approximation, the characteristic velocities of channels and ridges in the 
spacetime diagram of the probability distribution, and in section HVI we examine the way in which the degeneracy of 
these velocities shapes the quantum carpet. These ideas are illustrated with the example of the infinite square well 
potential in section Ivl and are used to elucidate the structure of fractional revivals in section IVTl Section IVlII offers a 
summary of lessons learned. 



II. MULTIMODE INTERFERENCE AND CHARACTERISTIC VELOCITIES 



A multimode term is a straightforward adaptation of a density matrix element, 

dnm = CnCjyi = (5) 

which allows us to write a probability density as 

oo 

= J2 f^nm{x,t). (6) 

n,m— 1 

While the multimode terms lack some of the convenient properties of the density matrix they can simplify the study 
of some dynamic phenomena, since they allow us to write the probability density as a simple sum. Of course, in order 
to do much with them we will have to make an assumption about the eigenfunctions. 

We will use the WKB approximation jl^] to obtain a usable form of the eigenfunctions, ipn- While this limits the 
range of applicability of our findings, quantum probability distributions only acquire the features that earn the name 
"quantum carpet" in the semiclassical case j30{ . The approximate eigenfunctions are 

^(— ) / r'^ \ /^(+) 



ipn{x) w , exp i / pn {x') dx' + exp ( ~i / p„ {x') dx' ] , (7) 

^yPn(X) \ J J \/Pn{x) 



where the ci^^ are complex constants of integration, and 



Pn{x) = ^En-Vix). (8) 

The approximate multimode terms fj,nm will be weighted sums of four intermode terms, 

a Kid 1'"")^ (9) 

where 

4™^'^(a^, t) = dn™-^=i^= exp ~i( [ (±ip„ (x') ±2 p^ {x')) dx' + {E,, ~ E„,) t] . (10) 

\/Pn\x)pm{x) \J J 

We are interested in features of the probability density that are time-dependent and typically much more "narrow" 

— 1/2 

than the potential, so we focus our attention on the phase factor instead of the leading factors of p„ ' . We can 
characterize the curves of constant phase by differentiating the argument of the exponential and looking for its roots, 

{±lPn{x')±2Pm{x'))dx' + [En- E^)t\ = 0, (11) 

r/r E — E 

J^r7, (12) 




±lPn [x) ±2 Pm (x) 



Afc ±,^E~V{x) ±2 ^E~V(x)' 

we find that each pair of quantum numbers, (n, m), gives rise to four velocities wim^^^^ We can pick any one of these 
velocities, integrate it from an arbitrary xq, and the resulting trajectory is a line of constant phase. Channels and 
ridges in the probability distribution follow these lines, as we will demonstrate in section Hvl 
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III. CHARACTERIZATION OF THE VELOCITIES 

We have written eq. ^] as Aoj/Ak anticipating that we wiU label some of the velocities as "group velocities." 
Defining a;„ — i?„ and fc„ (x) = p„ (x) , we can write 

LJn{kn,x)^kl-V{x). (14) 

We now make the additional assumption that the weighting coefficients are well-centered on some number n with 
some spread An, and that these satisfy the hierarchy 1 ^ An ^ n. This allows us to define a group velocity for our 
packet, 



duJn 



= 2fcfi, (15) 



which is the classical velocity of a particle in the potential V with energy En ■ 

Writing En = En + e„ and E^ = En + em, where e„,em ^ En — V(x), we can simplify our expression for the 
velocities in eq. El 



,.(±1±2) ^ _ 

nm 



[En + e„) - [En + e™) 



±i^jEn'Ven~V{x) ±2 yjEn + 6™ - V(a;) 

V^^-n^^) (±1 (1 + 2iE:-Vi.)) ) ±2 (1 + liE^^)) 

-2^En-V{x)- 



^n ^m 



2 {En - T^(a;)) (±il ±2 1) + (±ie„ ±2 e„) 



f ±2/E~y(^ = ±i2fcfi(a;), (+1-2, -1+2) 

^ 1 ^S^^^n^^^'^W' (+1+2,-1-2). (1 ^ 

As long as both e„,em ^ — ^(2;), s-s is required in the semiclassical approximation, then half of the velocities 
contributed by a particular (n,m) are comparable to the group velocity, Vgr, and half are smaller. It is important 
to note (in preparation for the next section) that eq. 1161 is approximate - it is not the case that half of the velocity 
terms are exactly degenerate. 

IV. GROUPS OF VELOCITIES AND DEGENERACY 

The result in eq. El demonstrates that in any problem there exists a range of velocities. Though that result was 
approximate, we will now show that some intermode terms have exactly the same maximum velocities, and that 
because of this we can treat those intermode terms as essentially moving together. This is a sort of degeneracy - the 
sort of degeneracy that produces quantum carpets. 

If we examine the velocities (eq. I13II at a point where V{x) — 0, we find that 

„(±i±.) = (17) 

±lV-C/n ±2 V £/m 



We can factor the numerator, 

^En + V Em ) ( \/ En — vE^ 



V. 



(±1±2) ^ 

nm I / 771 _|_ / 77' 



(18) 



and arrive at an important condition, 

vit^''^ ^Ti\/K±2\/e^- (19) 
Terms that have the same Vnrn at V{x) ~ change phase with the same period - these are our degenerate terms. 



We define a group of velocities as 



p.{x,t)= V 4±;±^)(x,t). (20) 



„(±1±2)| 
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Because the first sign in the definition of the velocities (eq. I19II is flipped relative to the signs in the definition of an 
intermode term (eq. I1(J|) . each /3^will look roughly like a weighted sum of harmonic functions. If there are many terms 
in the sum, they may interfere to produce a peak or a valley. As the wavefunction evolves in time, these peaks and 
valleys will trace out the canals and ridges that form the quantum carpet. 

This immediately demonstrates the role of quadratic spectra in producing quantum carpets. If the spectrum depends 
on the quantum number n squared, then there will be many degenerate velocities, while if the spectrum is linear in 
the quantum number, only traces which involve two perfect squares may be degenerate. This suggests that if we begin 
with a system like the simple harmonic oscillator, with a spectrum linear in the quantum number, and set to zero 
all weighting coefRcients that are not perfect squares (selecting only states 1,4,9,16, etc.), we can produce a carpet 
because we have effectively "quadratized" the spectrum [sj]. For an example of this, see Figure ^ 



V. THE INFINITE SQUARE WELL 



We now turn our attention to a particular example in order to illustrate the utility of these results. We write 
explicitly all of the relevant physical constants, so as to facilitate comparison with classical results. For a particle of 
mass M in a well of width L we have the wavefunction. 




We can immediately identify the velocities in question as 

= (=Fi" ±2 m) = vo (Tm ±2 m) . (22) 

We can also find that the wavenumbers knm will be 

fc±/„±^ -^(±in±2m), (23) 

so that in any particular velocity bundle there will be a variety of wavenumbers. We provide an example of the 
interference that this produces in figure |21 
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Figure 2: The |«„m| — vq (left) and |u„mj ~ 2uo (right) velocity bundles for a uniform distribution of weighting coefficients 
between n = 1 and n = 10 in the infinite square well. 




0.2 0.4 0.6 0.8 1 0.2 0.4 0.6 0.8 1 0.2 0.4 0.6 0.8 1 



X Lfi X Lfi X Lfi 

Figure 3: At left, the full carpet, coefficients evenly weighted between n = 1 and n — 10. In the middle, the eighteen intermode 
terms that compose /3„„. Note that they do, indeed, have a period T = Tr. At right the 34 intermode terms that compose 
l32vQ- Note that they have a period T = Tr/2. 



The most degenerate velocity will typically be vq, and in time T/f = 2L/vq a trajectory with this velocity will cover 
a distance of 2L, so the most prominent traces should have the same period as the revival time. Better still, we can 
ask what velocity we need in order to have a period equal to the classical period, T^. The condition is 

nVQ 

Vnra = 2riwo, (24) 

which should be satisfied by very few (n, ni) pairs. An example of how well our separation of the wavefunction works 
is shown in Figure |21 

In the case of the square well we can find a closed form for Pv{x, 0). From the definition and the form of uini it 
can be shown that 

Coo * oo ^ 

n— 1 m—l 

n=l / 

and recalling that the ti^^^'' are just exponentials, we can write (3y{x, 0) as 

2 °° TTX 1 'KX 

Pvix,t) = -^Re{c„<+„}cos(^(2n + u)^) +I'^°^(^t) E'^"'^^-"' (^^^ 

n—1 n—1 
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The first term contributes peaks and valleys, while the second sets a threshold height above which those peaks and 
valleys must rise in order to be visible. The first term suggests that we can, by choosing appropriate c„, control the 
distribution of peaks and valleys, though we cannot generate an arbitrary function [s^ . 



VI. MULTIMODE TRACES AND QUANTUM REVIVALS 

In problems where the weighting coefficients c„ are well-localized around some central quantum number n, so that 
the spectrum can be well approximated by 

1 2tt 2tt 

En^En + E'^ {n-n) + -E'^ {n - rl) = E^ + —k+ —k\ (26) 

^ J^cl -iR 

we can make an interesting approximation to the wavefunction at rational fractions of Tr. We do not present the 
details of the derivation here (see or [l^), but quote the result that 

^ (x^t^^TRj (2^' * + 7^-0 ' (27) 

where I — q/2 when q contains more than one power of two, I = q otherwise, the as are complex weighting coefficients 
(for their values, see or [2^), and we have defined 



fc— — 00 



^cz ix,t) = Ck+n^k+n{x) exp ( -2mk^^ ) . (28) 



This approximation states that the wavefunction can be written as the sum of slices of a "classicized" wavefunction, 
^ci- It suggests that at rational fractions of (particularly fractions with low g), a wavefunction that began as a 
well- localized packet will consist of several packets. We will now demonstrate that if \E'(x, t = 0) is well-localized, then 
^'c/ will remain well-localized for all t. 

Observe from the derivation of Equation El that the En ~ Em term comes from the time-evolution exponential and 
the square root terms come from the WKB approximation. If we want to find a similar formula for '^cU we need 
only replace the En — Em in the numerator with (27r/Ti) [n — to). Considering a point where V{x) — 0, we find the 
following velocity degeneracy condition: 

27r n — m 

w„„ = ±1 — . (29) 

Unlike Eauation ll9l we have no hope of factoring this in general. We can, however, see what happens in a few obvious 
cases. If 

En = o?n, (30) 

then we quickly find 

27r n-TO 

t;„™ = ±i- ^- — -j=, (31) 

icia yjn ±2 VTO 

a degeneracy condition identical to the one we would have found for the harmonic oscillator in the previous section. 
If, however, we try a spectrum En ~ o?r? ^ we find something far more interesting, 

_ 27r n-m _ \ ±1^, (^2) , , . 

Te;an±2TO I ±lT,,an-(-m' (. + 2) ■ 

Half of all of the traces are degenerate! Looking at the other half, we want to find two pairs of numbers, (n, m) and 
(p, q) that will be degenerate. The conditions are 

n-m _ p - q n + m=^0, p + q=^0, , . 

n + m^p + q'n^p, m ^ p, 







Tel 
2 
tfi 



cl 
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(a) The Morse Potential 



(b) The Rosen-Morse I Potential 



Figure 4: Plots of 'i'd in two potentials, one with a quadratic spectrum and one with a non-quadratic spectrum. On the left is 
the Morse Potential, with En = — {A — anf' , and on the right the Rosen- Morse I potential, with En = —A^ + {A-\- an)^ — 
{B/{A + an)y + {B/Ay. In both cases A, B, and a are independent constants. The black lines overlaying the maxima are 
classical paths for those potentials with E — En- These are exactly solvable potentials taken from [lo|. 



which then reduce to 



{p + l) (" — m) = {n + m) {p — q) , 

(34) 



qn — mp = —qn + mp, 
n p 

m q 



While this is not a difficult equation to satisfy with all of the integers at our disposal, it is difficult to satisfy when 
our weighting coefficients are all within some An <g; n of n 3> 1 . The solutions that will lie closest to (n, m) are 
(p, q) = 1/2 (n, m) and {p, q) = 2 (n, m). If n and m are in the vicinity of n, as they must be in the semiclassical 
case, that puts p and q near either n/2 or 2n, both of which would typically be "out of range" of An. For the case 
of a quadratic spectrum, then, we find that half of the velocities are degenerate, and that the other half tend to be 
non-degenerate. This nicely corresponds with our suspicion that wavepackets evolve along classical paths. 

The unfortunate aspect of this result for ^'c/ is that it is so heavily dependent on the quadratic character of the 
spectrum. This behavior seems to occur in potentials with non-quadratic spectra as well, as demonstrated in Figure 
IP However, in the semiclassical limit (that is, for very large n), their spectra can appear to be no more than quadratic 
|24j . This leaves low-n cases to consider, and completing our understanding of 'i'd is a possible subject for future 
work. 



VII. SUMMARY 



We have seen how an analysis of the intermode terms in the probability density of a semiclassical wavefunction can 
give us insight into the degeneracy conditions which must be met in order to produce a carpet, and into the carpet 
itself. When those degeneracy conditions are met, we can examine the initial form of certain bundles of velocities 
and from that find the channels and ridges of the quantum carpet. For a problem with a spectrum linear in the 
quantum number it was almost impossible to generate a carpet, though a perverse choice of weighting coefficients 
could effectively quadratize the spectrum. Similar reasoning suggested why '^ci should so often resemble a classically 
oscillating packet, and thus why fractional revivals of a wavepacket should themselves resemble the initial packet. 
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Appendix A: DEMONSTRATION THAT Ti < r2 < T3... 



We begin by assuming that our spectrum can be written as a polynomial in n, 

M 
m— 1 

In this case, the jth derivative of £"„ is 

and if n > TO — j the jth derivative will reduce the contribution from the ?ith term. So long as n 3> M, it will be true 

that E^n'^ > > E^n\-; since < 10^ < 1- Consequently, Ti < T2 < Tg... holds. Our resuh is not significantly 
changed by allowing negative powers of n. Consider 

M 

En = a™"""- (A3) 

m— 1 

Its derivatives are 

m=l ^ ^ ' ^ 

SO we still get that Ti ^ T2 ^ T3... so long as n ^ M. Transcendental spectra still present a problem, though we 
note that the spectrum of any one-dimensional potential cannot increase faster than rt} for large n, roughly because 
the walls of the potential cannot be "harder" than the infinite square well ■ 
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